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UNDAMPED  WAVES  IN  A  COLLISION  FREE  ELECTRON 
PLASMA 

0  0  0 

By  Donald  Me  |-EAN 


INTRODUCTION 


Recently  Bernstein,  Green  and  Kruskal  £"l  J  Have  shown  by  a  non-linear 
theory  that  there  can  exist  undamped  longitudinal  waves  in  a  "collision  free" 
plasma,  and  further,  they  claim  that  these  waves  do  not  satisfy  a  dispersion 
relation. 

This  is  in  complete  contradiction  to  the  results  of  Landau  i  2]  and  of 

Van  Kampen  £3jf  ,  both  based  on  the  same  linear  approximation  to  the  Boltzmann 

equation,  and  apparently  equivalent,  which  predict  that  all  sufficiently  smooth 

initial  distribution  functions  give  rise  to  damped  potential  waves,  the  minims* 

damping  rate  being  given  by  Landau.  It  is  true  that,  if  the  distr  bution  functions 

giving  rise  to  the  B.G.K,,  undamped  waves  ire  replaced  by  singular  distributions, 

in  the  linear  theory 

one  can  demonstrate  undamped  waves/  but  since  the  distribution  functions  corres¬ 
ponding  to  the  B.G.K.  equilibria  are  not  in  general  singular,  the  linear  and 
non-linear  theories  must  be  considered  to  be  in  contradiction. 

Landau  and  Van  Kampen  did  find  a  dispersion  relation  in  the  sense  that 
the  most  slowly  damped  waves  are  those  which  satisfy  the  relation,  previously  given 
by  Vlasov  (j4J. 

This  dispersion  relation  is  basically  the  same  as  that  found  by  Bohm 
and  Gross  [5j ,  by  assuming  that  there  are  no  particles  near  the  wave  velocity. 

In  the  same  paper  Bohm  and  Gross  showed  that  the  presence  of  trapped  electrons 
would  modify  their  results,  but  argued,  intuitively,  that  only  those  waves  satis¬ 
fy*'^  their  dispersion  relation  would  be  strongly  excited. 

Recently  Denisse  (jsj  ,  in  a  paper  largely  aimed  at  criticizing  Landau's 
results,  has  derived  an  expression  for  k,  the  wave  number,  which  must  be  satisfied 
by  small  amplitude,  stable  waves  and  which  is  in  form  identical  with  the  dispersion 
relation  of  Vlasov,  but  which  is  to  be  interpreted  slightly  differently,  due  to  its 
different  derivation.  Denisse 's  approach  is  basically  similar  to  that  of  Bernstein, 
Greene  and  Kruskal  and  his  linear  approximation  is  more  satisfactory  than  that  of 
the  linearised  Boltzmann  equation  in  that  he  only  neglects  small  quantities. 

Most  of  this  report  is  concerned  with  an  examination  of  the  properties 
of  such  small  amplitude,  stable  waves  based  on  Denisse 's  results,  which  are  quoted 
in  Section  III. 

Before  this  in  Section  II,  a  short  discussion  is  given  of  the  bearing 
of  the  paper  of  Bernstein,  Greene  and  Kreskal  on  the  theory  of  Landau.  In  parti- 


cular  an  exanple  of  a  wave,  for  which  the  calculations  can  be  completed  exactly, 
and  which  is  undamped,  is  presented. 
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In  Section  III,  after  the  presentation  of  Denisse's  results  the  intets. 
pretation  and  significance  is  c  scussed. 

In  Section  IV  this  discussion  is  continued  to  demonstrate  the  absence 
of  a  true  dispersion  relation,  although  for  a  different  reason  from  that  given 
by  Bernstein,  Greene  and  Kruskal. 

In  Section  V  the  results  of  section  IV  are  used  to  criticize  the  results 

of  Landau. 

Finally  in  Section  VI  an  attempt  is  made  to  consider  the  effect  of  colli¬ 
sions  within  the  framework  of  the  theory  of  stable  solutions.  Such  an  approach  is 
of  course,  impossible  to  follow  up  rigorously,  but  the  somewhat  intuitive  results 
seem  sufficiently  interesting  to  be  worth  presenting,  as  thay  lead  to  a  new  inter¬ 
pretation  of  the  dispersion  relation. 
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SECTION  II 

Very  briefly,  Bernstein,  Greene  and  Kruskal* s  method  consisted  of  assu¬ 
ming  a  stationary  wave  and  finding  a  differential  equation  which  the  potential,  (p  , 
must  satisfy  this  assumption  (collisions  are  being  neglected).  This  equation  can 
equally  well  be  interpreted  as  an  integral  equation  for  the  distribution  function 
of  the  trapped  electrons  as  a  function  of  energy  -if <p (x)  and  the  distribution 
functions  for  protons  and  untrapped  electrons  are  specified.  Bernstein,  Greene  and 
Kruskal  have  presented  the  solution  of  the  equation  in  this  form. 

It  is  not  however  completely  clear  that  these  solutions,  the  existence  of  which 
is  the  essential  point  of  their  argument,  can  be  physically  meaningful.  Montgomery 
and  Gorman  for  example,  question  the  analyticity  of  these  solutions.  It  is 
therefore  interesting  to  demonstrate  a  particular  case  of  a  stationary  wave  for 
which  the  distribution  function  is  an  analytic  function  of  the  position  x  and  the 
velocity  v  (for  x  and  v  real).  In  the  Appendix  A,  it  is  demonstrated  that  at  least 
one  such  solution  exists.  The  electron  distribution  function  considered  is  « 


f  (x  ,  v) 


_ V _ 

v2  +  v2  +  2a^(x) 


where  n  is  the  electron  density  at  the  potential  maximum,  at  which  point  the  po- 
oe 

tential  <|>  (x)  is  zero,  V  is  an  arbitrary  positive  constant,  q^  and  m  are  respecti¬ 
vely  the  electronic  charge  and  mass  and  the  potential  (x)  is  an  analytic  func¬ 
tion  of  x,  the  inverse  of  which  is  derived  analytically  in  the  appendix. 

As  will  be  seen  in  the  Appendia  A  this  is  not  a  particularly  interesting 
example  physically  t  not  only  is  the  temperature  of  the  plasma  infinite,  but  the 
potential  (x)  tends  to  -  Oo  exponentially  when  x  tends  to  £  OO  •  Nevertheless 
it  proves  that  not  all  the  B.G.K.  equilibria  have  singular  distribution  functions 
(as  a  function  of  x  and  v). 

It  has  been  shown  by  Bernstein,  Greene  and  Kruskal  that  their  non-linear 
solutions  correspond  to  singular  approximate  distribution  functions  in  the  linear 
approximation  but  this  does  not  "salvage  the  conventional  theory"  since  this 
implies  that  before  performing  the  linear  calculations  we  must  replace  the  true 
initial  distribution  function  by  another  modified  distribution,  but  we  only  know 
how  to  do  this  (or  even  that  it  is  possible)  for  the  case  of  the  B.G.K.  equilibria. 
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The  paper  of  Bernstein,  Greene  and  Kruskal  does  not  attempt  to  treat  the 

geperal  problem,  treated  by  Landau,  of  describing  the  evolution  in  time  of  a  plasma 

assumed 

which  is  perturbed  at  the  xime  t  ■  0.  More  exactly  Landau/the  initial  distribution 
f  (x,  v,  0)  known  and  sought  to  derive  the  distribution  ,f  (x,  v,  t)  subject 
to  the  condition  that  f(x,  v,  0)  represents  a  small  perturbation  from  the 
equilibrium  conditions.  He  found  that  the  electric  field  E(x,  t)—^  0  when 
Effectively  Bernstein  et  al.  have  produced  a  counter  example.  If  we  take  the  special 
case  in  which  f{x,  v,  0)  is  the  distribution  corresponding  to  an  undamped  wave  we 
must  have  f(x,  v,  t)  »  f  (x  -  ut,  v,  0).  Whence  E(x,  t)  =>  E(x  -ut,  0)  /-^  0. 

(u  is  the  velocity  of  the  wave  relative  to  the  frame  of  reference). 

Hence,  subject  to  the  assumption  that  the  equations  used  have  a  unique 
solution  (which  has  apparently  been  proven  by  Iordanskii  CeJ  ).  Landau's  results 
are  disproven  by  those  of  Bernstein  et  al. 

The  rest  of  this  report  is  concerned  with  the  properties  of  stable  plasma 
waves  -  the  existence  of  which  is  not  further  questioned,  except  for  a  short  demons¬ 
tration  in  Section  V  that  Landau's  results  depend  critically  on  hi3  assumptions. 
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Those  of  Denisse’s  results  which  interei t  us  especially  here  can  be 
collected  into  the  following  three  equations  t 


q.  +  n  q 
1  oe  e 


+  o  ( 


i  '“o’ 


r  ,v  v 
o  (  o  ) 


Here,  x  is  the  distance  coordinate  in  the  frame  of  reference  in  which 

the  wave  is  stationary,  parallel  to  the  direction  of  propagation, Cp (x)  is  the 

electric  potential,  n.  the  ion  density  (assumed  constant)  q,  and  q  the  ionic  and 

i  l  6 

electronic  charge  respectively,  m  the  mass  of  an  electron.  The  potential  is  assumed 

to  have  a  maximum  which  can  be  taken  os  zero;  v  is  then  the  velocity  of  an  electron 

o 

along  the  x-axis  at  this  maximum  and  f  (v  )  the  electron  distribution  function 

point  °  o 

at  this  point  (which  is  the  only/through  which  all  electrons  pass). 

(k^  is  the  permittivity  of  free  space). 

These  equations  have  been  derived  on  the  assumption  that  all  quantities 
are  independent  of  time  in  some  frame  of  ref  ren'-e,  and  have  no  sense  otherwise. 

Further,  it  must  be  assumed  that  ~  n-  is  continuous  for  v  =0, 

avo 

otherwise  the  Cauchy  principal  value  of  the  integral  indicated  in  equation  (2)  does 
not  exist.  This  does  not  appear  to  be  an  unreasonable  assumption  physically. 

We  can  at  once  write  the  solution  of  equation  (1),  valid  for  ^mgx 

where  i8  some  arbitrary  level  above  which  the  higher  order  terms  in^cannot 

be  neglected. 

We  have  for  fc2  0 

M  -  j^"cos  k(x  -  xQ  )  -  jTj  (4) 


and  for 


K2  <  0 


<J>(x)  «  -  |cos^K  (x  -x0)  - 
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where 


n  q  +  n  q 
oe  "e  i  i 


(6) 


The  assumption  that  ^23  requires  that  °^/0  ^  0 

ie  n ,  \  n 

i  oe 


Clearly,  in  equation  (4),  n.  can  be  chosen  sufficiently  close  to  n 

A  1  oe 

that  — »- 

and  in  this  case  equation  (4)  will  be  a  good  approximation  to  the  potential  which 


Q)  (the  limit  of  validity  for  the  linearized  form  of  equation  (1), 
y  max 


could  be  predicted  by  Bernstein,  Greene  and  Kruskal's  non-linear  theory,  if  it  were 

2 

calculable.  Hence  if  k  f  derived  from  equation  (2)  is  positive,  it  yields  the  wave 
number,  for  small  (potential)  amplitude  waves.  On  the  other  hand,  since  cosh  x  i6 
unbounded  above,  equation  (5)  can  only  ever  be  a  good  approximation  for  a  small 
range  of  x.  In  fact,  in  this  case,  the  potential  is  not  even  necessarily  periodic 
and  so  knowledge  of  the  value  of  K  is  of  very  little  value. 

When  the  right  hand  side  of  equation  (2)  is  negative  there  are  no  waves 
of  arbitrary  small  amplitude  -  but  beyond  this  equation  (2)  cannot  be  considered 
as  giving  any  useful  information.  In  particular  it  should  not  be  interpreted  as 
indicating  evanescent  waves,  such  waves  being  explicitely  excluded  from  Denisse'3 
calculations. 

On  the  other  hand,  when  the  right  hand  side  of  equation  (2)  is  positive, 

equation  (2)  has  the  sense  of  dispersion  relation,  except  that  we  have  not  yet 

specified  f  »  Such  a  specification  requires  a  detailed  description  the  menntfs  in 

which  the  wave  is  excited  .  Such  a  description  is  of  course  beyond  the  scope  of 

any  theory  so  far  developed.  The  only  logical  procedure  is  to  consider  fQ  as  an 

arbitrary  function,  the  different  functions  f  corresponding  to  different  mecha- 

o 

nisms  of  excitation.  This  question  is  taken  up  in  detail  in  the  next  section. 

It  should  be  noted  however,  that  physically  we  must  exclude  distribution 
functions,  f  ,  which  are  not  symmetric  about  zero  for  all  velocities  less  than  that 
nec*»«3ary  for  an  electron  to  "escape"  from  the  potential  trough.  In  the  limit  of 
small  amplitudes,  in  which  we  are  interested  here,  this  interval  of  velocities 
vanishes,  and  it  is  not  unreasonable  to  ignore  the  restriction,  provided  we  retain 
the  symbol  for  the  principal  value  of  the  integral  in  equation  (2).  Even  if  we  were 
to  extend  the  theory  to  waves  of  non-zero  amplitude,  we  can  formally  avoid  this 
restriction  because  the  contribution  of  an  electron  to  the  charge  dkensity  (through 
the  time  it  spends  in  each  part  of  the  wave)  depends  only  on  its  speed  relative  to  the 
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wave |  and  not  on  the  direction  of  its  velocity,  we  find  that  there  is  no  change  to 

the  wave  if  we  reverse  the  velocity  of  a  group  of  particles.  This  is  true,  both  in 

the  non-linear  theory  of  Bernstein  e^jal.  and  in  Denisse's  linear  approximation. 

Consequently  (as  can  be  easily  verified  by  substitution  in  equation  (2)  ,  the  results 

obtained  for  a  given  f  (v^)  are  identical  with  those  obtained  for  any  distribution 

which  can  be  constructed  by  adding  an  odd  function  to  f  (v  )  -  including  of  course 

oo 

a  set  of  functions  which  are  symmetric  in  any  required  range.  Hende,  when,  in  what 

follows,  we  calculate  results  for  functions  f  (v  )  which  do  not  satisfy  the  sym- 

00  all 

metry  conditions  in  any  range,  we  are  really  calculating,  fat  those  functions, 
symmetric  in  the  required  range,  which  could  be  constructed  by  the  addition  of  an 
odd  function.  In  the  limit  of  small  amplitudes,  in  which  we  are  interested  here  the 
difference  is  very  small. 
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SECTION  IV 

In  this  section  wc  discuss  what  significance  can  be  given  to  equation  (2) 
in  the  light  of  the  point,  made  in  the  last  section,  that  f  (v^)  must  bo  considered 
as  an  arbitrary  function,  unless  the  excitation  mechanism  can  be  described. 

Bernstein,  Greene  and  Kruskal  argue  that  no  dispersion  relation  can  exist, 
due  to  the  fact  that  the  velocity  (relative  to  the  wave)  of  groups  of  untrapped 
electrons  can  be  reversed  without  modifying  the  potential.  This  has  already  been 
pointed  out  in  the  last  section  and  is  represented  in  their  calculations  by  the  fact 
that  it  is  only  necessary  to  know  the  distribution  of  the  particle  energies  (measured 
in  the  wave  frame)  to  calculate  the  potential.  Using  this  fact  the  velocity  of  the 
plasma  relative  to  the  wave  or  the  wave  velocity  relative  to  the  plasma,  can  be 
modified  by  reversing  the  velocity  relative  to  the  wave  of  some  of  the  electrons 
without  in  any  way  modifying  the  potential  wave,  and  so  there  is  no  relation  between 
wave  velocity  and  wavelength,  or  between  frequency  and  wave  number. 

However,  in  seeking  a  dispersion  relation,  one  expects  such  a  relation 
to  reflect  the  nature  of  the  unperturbed  plasma  -  for  instance  Bohm  and  Gross  derived 
on  approximate  dispersion  relation  involving  the  density  and  temperature  of  the 
plasma.  If  however  the  disturbance  of  the  plasma  is  very  great,  the  temperature, 
for  instance,  may  be  greatly  modified,  in  which  case  the  dispersion  relation  calcu¬ 
lated  using  the  "equilibrium"  temperature  is  no  longer  obeyed.  Hence  we  can  only 
expect  a  dispersion  relation  to  have  a  sense  if  the  deviation  from  equilibrium  is 
small,  and  it  is  not  surprising  that  Bernstein  et  al.  did  not  find  a  dispersion 
relation  without  imposing  such  a  restriction. 

In  this  section  we  consider  the  existence  of  a  dispersion  relation  for 
stable  waves  subject  to  the  additional  condition  (which  clearly  limits  the  argument 
used  by  Bernstein  et  al.  thsti 


f  (x,  v) 


f  (v  -  v  ) 
eq  p 


(7) 


is  satisfied  for  all  x  and  v  I  v^  is  effectively  the  phase  velocity  of  the  wave,  and 


is  considered  as  a  parameter,  f  is  the  Maxwell  distribution  for  the  velocity 

eq 

component  parallel  to  the  direction  of  propagation  of  the  wave  defined  by  equation  (0)  i 

/ - -n  _  —  a--  V2  (0) 

f  (v)  -  n  (-  ...7! -  e  2kT 

«  s/  2  7T  kT 

(  k  <  Boltzmann's  constant;  Ti  electron  temperature;  n^  ;  equilibrium  electron  density). 
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From  equation  (7)  we  deduce  that  the  amplitude  of  the  potential  must  be 
small,  since  otherwise  the  density  varies  greatly  from  its  equilibrium  value  and  (7) 
is  certainly  not  satisfied.  Hence  we  may  apply  equation  (2).  In  addition  we  shall 
make  use  of  a  special  case  of  equation  (7)  s 


f  (v  )  f  (v  -v  ) 

o  0  eq  o  p 


(7») 


which  we  rewrite  f  (v  ) 
o  o. 


D(f  (v  -  v  )  +  Pif(v  ) 
eq  o  p  o 


(7") 


It  should  be  noted  that  f  (v)  is  the  equilibrium  distribution  in  the 

eq  • 

rest  frame  of  the  plasma  but  f  (v  )  is  the  distribution  at  a  point  in  a  particular 

o  o 

wave,  the  velocities  being  measured  in  the  frame  in  which  the  wave  is  stationary. 

Here# is  approximately  unity  and  £}f(vQ)  is  small  compared  with 

f  (v  -  v  ). 
eq  o  p 

It  is  possible  to  examine  a  number  of  special  cases  analytically,  in  order 
to  study  the  behaviour  of  equation  (2)  subject  to  the  conditions  of  equation  (71). 

The  first  i«  the  obvious  case  ^f  5  0. 


f  (v  )  *  Ct  f  (v  -  v  ) 

o  o  wv*  eq  o  p 


(9) 


Combining  equations  (2),  (8)  and  (9),  we  obtain  s 

,  00  2  1 
,2  P  1 


where 


H2 


2  Vtf 

0C\%  "Bq, 


e  e  ^  ene 


K  m 
o 


K  m 
o 


the  equilibrium  distribution, 


kT 


91 


P  ) 


(10) 


the  square  of  the  plasma  frequency  of 


-oo 


and 


g,  (y) 


x  e 


x  -  y 


dx 


(ID 


—  oo 


In  Appendix  B  it  is  shown  that  the  function  g^  (y),  defined  by  equation 
(11),  is  given  by  t 


(y) 


2  y  e 


-  y2 


J 


du  -  1 


(12) 
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This  function  is  tabulated  is  various  places  (see  reference  8  ).  A 

rough  sketch  of  k2  given  by  equation  (10),  deduced  from  various  approximations  which 
can  readily  be  obtained  from  equation  (12)  is  shown  in  Fig,  1, 

In  Appendix  C  it  is  shown  that  g^(y)  can  be  expanded  in  an  esymptotic 
'series,  of  which  the  first  two  terms  are  t 


h  (y)  fjf  (  “V  +  l 

2y 


(13) 


Hence 


.2  ,2  2  ?  . 
<*>  "  k  VP  SJp.iV2 


H2 

wT- 


OJp2  ..  ,  . 

In  the  region  for  which  this  expansion  is  useful  and  S°  WB  '*n“  'or 

the  case  of  small  k 


co2 


£2  ^p2  + 


3  kT 
m 


(13') 


which  is  the  Bohm  and  Gross  dispersion  relation. 
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The  approximation  to  k  ,  obtained  by  substituting  (13)  in  (10)  is 
presented  as  a  broken  curve  of  figure  1. 

We  now  consider  the  second  special  case 

2 


0f(v  ) 


6" 


\Z27V1  W 


Vn  _  W 

o'  1/2  (  “V- 


(14) 


The  resulting  distribution  function  f  (v  )  is  shown  in  figure  2. 

2  °  °  2 

The  correction  0(k  )f  which  must  be  made  to  k  which  we  have  already 
deduced  from  equation  (10)  is  simply  : 


8(k2) 


i±3aS. 


v9fw 


_i_ 

w 


/  —  VXr- 

1  #w 


) 


(15) 


This  correction  is  plotted  in  figure  3  os  a  function  of  w  for  several  values  of  W 
It  is  of  course  the  same  curve  as  figure  2  on  a  different  scale,  with  different 
coordinates. 

An  examination  of  this  curve  brings  out  the  following  points.  If  W  is 

large  0(k2)  is  small  for  all.  w.  If  w  is  large,  0(k^)  is  small,  that  is  to  say 

small  modifications  to  f  (v  )  do  not  affect  k  in  the  region  where  v0  0. 

o  o 

More  precisely  t 
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&  (k2)  k2 


if 


LQ. 


2  GJ2 
S - a— 

w2- 


-  «  i 


(16) 


Since  this  .  inequality  doss  not  involve  the  form  of  0  f ,  it  is  not 
un  reasonable  to  suppose  that  it  holds  for  localized  modifications  to  f0  other 
than  gaussian  ones. 

On  the  other  hand  we  see  that  for  w  0, 
from  zero.  Using  the  properties  of  g^  (x)  (g-f  (0)  »  9^  (x)  ^  v/TT  ,wc  3eB 

from  equation  (15)  that 


0(k^)  may  differ  greatly 

■>  c 


mnx 


The  condition  (7f)  (i.e.fQ  *0^ 


6n.  ,  yg 
n«s  \/t?W  \V 

,  )  requires  that 


(17) 


be 


6n 

Jr  2  f 

£>n.  Hence/0(k  )/ 


small  but  this  condition  cc.n  be  satisfied  for  any  W  by  chodsing 
can  be  made  as  large  as  desired  by  reducing  W  in  equation  (17)  at  the  same  time 
reducing  0n  so  that  — ~- 


max 


W 


remSins  constant. 


We  see  at  once  that  the  condition  of  equation  (7)  (f  (x,  v);Njf  (v-v  )) 

oq  p 

is  not  sufficient  to  give  to  equation  (2)  the  meaning  of  a  true  dispersion  relation. 

We  also  see  another  important  point.  That  is  that  a  crude  knowledge  of 
the  distribution  of  those  electrons  which  move  very  rapidly  relative  to  the  wave  is 
sufficient  to  estimate  their  contribution  to  the  wave  -  but  that  the  wavelength 
(and  in  fact  the  form)  of  the  wave  is  very  sensitive  to  the  precise  details  of  the 
distribution  of  those  electrons  (trapped  or  almost  trapped)  which  move  slowly  rela¬ 
tive  to  the  wave. 

This  is  a  very  important  point  because  all  approximate  calculations  based 
on  tho  "linearized”  Boltzmann  equation,  as  well  as  the  approximate  calculations  of 
Bohm  and  Gross  -  give  an  accurate  description  of  those  particles  moving  rapidly 
relative  to  the  wave  (for  which  we  find  that  a  crude  description  suffices)  and  give 
a  quite  false  description  of  those  particles  with  velocities  close  to  the  wave  velocity 
(to  which,  as  we  have  just  seen,  the  wave  is  very  sensitive). 

It  is  possible  to  argue  that  thB  absence  of  a  dispersion  relation  can  be 
demonstrated  from  Bernstein,  Greene  and  Kruskal's  results,  since  given  the  distribution 
for  the  untrapped  electrons,  it  is  possible  to  derive  a  oist ribution  of  trapped 
electrons,  which  together  with  the  untrapped  electrons,  give  rise  to  any  desired 
potential  wave  (subject  to  certain  symmetry  conditions  when  we  ignore  the  iorfiotion) . 
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However  It  does  not  seem  possible  to  be  sure  that  the  derived  distribution  function 

will  not  be  very  great  or,  worse j  negative,  for  an  arbitrary  potential  wave.  It  therefore 

seems  best  to  argue  from  Denisse's  results* 

To  summarize  this  section,  if  we  substitute  ^eq^vo  “  vp)  f°r  fo^vo^ 

equation  (2)  we  find  effectively  the  same  dispersion  relation  as  Bohm  and  Gross  or 

(except  for  the  damping)  Landau.  It  does  not  however  seem  worth  while  discussing  the 

curve  representing  this  dispersion  relation  in  too  much  detail,  because  as  wo  have 

•  w 

shown  using  the  special  case  of  the  sum  of  several  gaussian  functions,  we  can  choose 

fQ(vQ)  to  differ  from  feq  <v0  “  vp)  by  as  little  as  we  please,  and  yet  k^  differs  from 

the  dispersion  relation  by  as  much  as  we  please.  In  other  words,  the  conditions  we 

have  imposed  on  f  (v  )  are  not  sufficient  to  ensure  the  existence  of  a  dispersion 
o  o 

relation.  The  question  remains  as  to  whether  physiaally  meaningful  additional  condi¬ 
tions  cannot  be  imposed  on  f  -  which  lead  to  a  dispersion  relation  in  some  sense. 

□ 

This  problem  is  taken  up  in  Section  VI. 


SECTION  V 
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In  this  section  the  value  of  a  theory  which  la  baaed  on  a  linear  approxi¬ 
mation  to  the  Boltzmann  equation  is  questions  1  (in  particular  Landau's  theory  and, 
aince  it..is  equivalent,  that  of  Van  Kampen). 

We  have  already  seen  in  Section  II  that  Bernstein,  Greene  and  Kruskal's 
non-linear  stationary  solutions  serve  as  counter  examples  to  Landau's  theory. 

As  pointed  out  in  Section  IV  for  the  case  of  stable  waves  the  linear 
Boltzmann  equation  describee  the  untrapped  particles  quite  well  but  the  trapped 
particles  very  badly.  In  fact  Landau's  method  amounts  to  treating  all  particles 
according  to  an  approximation  which  is  sound  for  those  moving  very  rapidly  relative 
to  the  wave.  On  the  contrary  Denisse  QfJ  has  been  able  to  make  the  distinction  bet¬ 
ween  these  two  groups  of  particles,  in  that  a  range  of  velocities,  corresponding  to 
some  of  the  trapped  particles,  are  always  excluded  from  his  integrals.  He  has  also 
shown  that  if  one  removes  this  distinction  between  trapped  and  untrapped  particles 
by  (quite  incorrectly)  integrating  over  all  velocities  the  extra  term  which  results 
in  the  "dispersion"  equation  has  exactly  the  same  form  as  Landau's  damping  term. 
Since  this  extra  term  corresponds  to  the  lack  of  distinction  between  trapped  and 
untrepped  electrons,  inherent  in  Landau's  approach  there  is  a  very  strong  suggestion 
that  it  is  from  this  aspect  of  his  approximation  that  the  erroneous  damping  term 
arises. 

In  Section  IV  we  have  seen  that  the  Landau  approximation  is  crude  where 
it  needs  to  be  precise  (trapped  particles)  and  precises  where  it  can  afford  to  be 
crude  (untrapped  particles)  which  in  itself  is  a  rather  severe  criticism. 

Finally,  in  the  rest  of  this  section  we  present  a  new  argument  against 
Landau's  theory  by  re-examining  his  assumptions.  Although  he  did  not  write  it  all 
down  as  explicitely,  his  argrment  starts  as  follows  (The  notation  is  necessarili 
different  from  Landau's). 

The  equations  to  be  solved  are  t 


O f(x,v,t) 
8  * 

-  *  V  ,  "»  E(..t, 

3  x  » 

3f(x,V,t)  n 

a  v 

da) 

0  E(x,t) 

a  x 

1 

r 

'  qe  J 

dv  f(x,v,t)  +  niqi 

(19) 

In  order  to  develop  a  perturbation  theory  one  writes  i 
f(x,v,t)  ■  Fq(v)  +  fj  (x,v,t) 


(20) 
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subject  to  the  assumptions  that  [  f^  | 


F  and 
o 


r  r 

aq 


which  implies  t 


F  ( v) 

O 


A7 


f  (v) 
eq 


(21) 


Landau  straight  away  took  F^(v) 


f  (v),  which  is  of  course  the  obvious 
eq 

choice,  but  not  the  only  one  valid.  We  shall  see  however  that  the  results  of  his 
calculations  depend  critically  on  the  method  in  which  the  division  between  the  zero 
and  first  order  distributions  is  carried  out  in  equation  (20).  It  is  of  course  just 
as  valid  to  carry  out  the  calculations  for  an  arbitrary  F^  subject  to  the  condition 
(21).  We  shall  however  add  here  the  condition  that 


F  (v)  dv  - 
o 


eq 


(v)  dv  =»  n. 


which  gives  for  equatim  (19) 

2Ll. 

ax 


'  Oo 


-T?” 

o 


f.  (x, v,t)  dv 


(19*) 


■oo 


This  is  convenient  in  that  the  remaining  calculations  are  identical  to  Landau's 
and  we  may  simply  take  over  his  results,  with  f  replaced  by  Fq  . 


Landau  then  " approximates"  to  equation  (IB)  by 

3'i 


3 1 


♦  *21  +  Jb. 

3x  "> 


E  ^ 

0  v 


the  term 


% 

m 


3  f. 


being  dropped  on  the  grounds  that  E  and 


3f1 


(IB*) 


are  both 


3*  —  —  *  3  v 

first  order  quantities.  Since  neither  the  integral  nor  the  derivative  of  a  small 

quantity  is  necessarily  small  this  approximation  is  immediately  questionable  and 

should  have  been  justified  a  posteriori.  In  fact  Landau  avoided  the  justification 

by  taking  the  equation  (10*)  und  (19*)  (with  F  ■  f  )  as  his  starting  point.  After 

o  eq 

this  starting  point  his  procedure  appears  to  be  valid,  as  is  born  out  by  Van  Kompen. 

We  can  therefore  use  his  results  to  demonstrate  that  his  implicit  assumptions,  stated 

explicitely  above,  are  unsound.  It  is  only  necessary  to  replace  f  by  Fq  in  all 

his  equations  and  note  the  effect  of  small  variations  of  F  on  G0(the  wavs  fre- 

o 

quency)  and  "Y^the  damping  rate). 

If  V  is  very  small  GO  i9  given  approximately  by  the  root  or  roots  of  s 
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and 


:Y  by 

Y 


'V, 


CO 


.7 T 


r  f*  (v) 


CO 


v  -  OJ/k 

2  F*  ( CO/k) 


dv 


dk 


(^2.) 

k 


(22) 


(23) 


where  k  is  the  wave  number  of  the  initial  perturbation. (See  for  example  Weenif*k 

For  a  derivation  of  these  results  without  assuming  Fj  «.  f  ) .  The  equation 

(23)  clearly  depends  on  the  assumed  value  of  F  ,  but  it  is  not  immediately  clear 

o 

what  the  effect  on  o  of  small  modification  to  F  will  be.  However  it  is  shown 

o 

in  Appendix  D,  that  we  can  choose  F  to  satisfy  the  following  condition  i 

o 


F  (v)  -  f  (v) 
o  eq 


<  e 


(24) 


for  any  £  0,  and  such  that  Landau’s  formulae  predict  waves  for  which  the 

damping  rate  is  identically  zero. 

Further  it  has  already  been  shown  in  Section  IV  that  equation  (22) 
does  not  relate  CO/k  (or  CO )  and  k  if  F^  is  only  subject  to  the  condition  (21). 
(Compare  (21)  and  (22)  with-  (2)  and  (7').) 

Landau's  results  therefore  depend  critically  on  his  assumption  that 
Fq(v)  m  f  f  but  this  assumption  is  quite  arbitrary. 

To  conclude  this  section,  we  see  that  the  following  arguments  can  be 
made  against  Landau's  paper. 

1 )  Bernstein,  Greene  and  Kruskal  have  produced  a  counter  example  which 
contradicts  Lg-idau.  This  in  itself  is  a  conclusive  point  which  continues  to  be 
ignored . 

2)  In  the  case  of  stable,  or  almost  stable,  longitudinal  waves,  the 
linear  Boltzmann  equation  gives  a  very  inaccurate  description  of  the  trapped  par¬ 
ticles,  which  are  just  the  particles  which  are  most  significant  in  determining  the 
wavelength  and  foim  of  the  waves. 

3)  In  Denisse's  theory  one  can  reproduce  Landau's  damping  term  by 
introducing  an  error  which  is  equivalent  to  the  error  inherent  in  a  linearized 
Boltzmann  equation. 
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4)  The  dispersion  relation  derived  by  Landau  (within  the  sense  he  gave  it) 
depends  very  critically  on  the  zero  order  distribution  function  on  which  the  solu¬ 
tion  is  based.  This  demonstrates  that  it  is  not  a  sound  perturbation  calculation, 
since  a  first  order  change  in  the  atro  order  function  should  not  have  an  enc.-nr.ous 
effect  on  the  results. 

Any  one  of  these  points  is  sufficient  to  raise  serious  doubts  as  to  the 
validity  of  the  results  of  Landau  and  of  all  equivalent  calculations. 
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SECTION  VI 

This  report  has  mostly  been  based  on  Denisse’s  results,  which  in  turn 
are  based  on  the  assumptions,  generally  made,  that  collisions  can  be  neglected, 
and  that  the  behaviour  of  a(large)  set  of  "point"  particles  can  be  adequately 
described  by  a  distribution  function. 

However  it  is  known  from  a  very  general  theorem  (Chapman  and  Cowling 
QoJ  )  that  in  the  presence  of  collisions  all  distributions  tend  towards  a  uniform 
Maxwellian  distribution  and  therefore  that  all  waves  are  damped.  It  is  generally 
assumed,  as  justification  for  neglecting  the  effect  of  collisions,  that  the  damping 
rate  due  to  collisions  will  be  of  the  order  of  one  of  the  relaxation  times  for  a 
"test  particle"  in  a  plasma  calculated  by  Spitzer  (for  example)  [~1l]  and  therefore 
negligible. 


In  this  section  it  is  planned  to  show  that  this  is  not  necessarily  so. 

It  is  of  course  impossible  to  study  rigourously,  within  the  framework  of  a  theory  of 
stationary  modes,  what  the  effect  of  collisions  will  be.  However  the  calculations 
which  will  be  presented  here  do  give  an  acceptable  physical  picture  of  the  signi¬ 
ficance  of  collisions. 

Before  proceeding  with  this,  however,  it  should  be  noted  that  the  repre¬ 
sentation  of  a  plasma  by  a  distribution  function  needs  careful  justification  -  and 
the  properties  of  a  distribution  function,  if  such  can  be  adequately  defined, 
should  be  carefully  examined. 

That  such  an  examination  is  necessary  can  be  seen  from  Van  Kampen’s^VJ 
classic  paper  on  the  subject  of  plasma  oscillations  in  the  absence  of  collisions. 

In  the  section  in  which  he  demonstrates  the  existence  of  normal  modes,  he  asserts 
that  a  distribution  function  is  sufficiently  well  defined  and  physically  accept¬ 
able  if  only  it  can  be  used  to  calculate  averages.  He  then  demonstretes  his  set  of 
modes,  which  are  singular.  However,  in  a  later  section,  when  he  decomposes  an 
"arbitrary"  initial  distribution  function  which  is  of  the  form  gQ(v)  e  into 
these  modes,  he  assumes  that  the  function 


9o+  (u)  ’  *  9o(u)  +  ~2~~l  P  - -  du' 

1  u  -  u’ 

does  not  include  any  delta  functions.  In  fact  this  excludes  singular  distributions, 
quits  contrary  to  the  argument  that  a  distribution  function  is  valid  as  long  as  it 
can  be  used  to  calculate  averages. 
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Turning  now  to  the  effect  of  collisions  we  consider  only  the  case  in 
which  the  time  of  interest  is  much  shorter  that  the  time  between  collisions, 
that  is  the  probability  of  a  given  electron  suffering  a  close  collision  and  we 
are  only  interested  in  the  effect  of  the  large  number  of  "distant  collisions", 
which  any  electron  suffers  is  even  a  very  short  interval  of  time. 

We  shall  further  limit  the  argument  by  considering  only  those  distribu¬ 
tion  functions,  of  the  form  considered  in  Section  IV  and  defined  by  equations  (7") 
and  (14),  i.e.  t 


f n  (V  ) 
O  o 


6f<“  > 


a-'., 

<0" 

Ww" 


-  V  +  Sf<v0> 

-  *  - )2 

e  W 


(7") 

(14) 


In  this  case  the  effect  of  collisions  on  the  function  f  will  be  ne- 

eq 

gligible  at  least  if  the  wave  amplitude  is  sufficiently  small  since  this  is  just 

the  function  which  is  not  affected  by  collisions  (in  a  uniform  plasma). 

We  are  left  with  the  problem  of  the  effect  of  collisions  on  £)f.  For 

this  it  is  convenient  to  replace  the  velocity  v  of  a  particle  by  the  velocity  vQ 

which  this  particle  would  have  at  its  next  passage  through  a  potential  maximum,  if 

there  were  no  collisions.  In  the  absence  of  collisions  v  woulu  be  a  constant  for 

o 

a  given  particle  by  definition.  In  fact  vq  will  vary  in  time  due  to  collisions,  in 

a  manner  which  it  is  impossible  to  determine.  It  is  however  possible  to  estimate 

the  probability  that  a  perticlp  which  has  a  velocity  v  at  time  t  will  have  a 

velocity  VQ  +  5  v0  a*  *  +&T.  In  fact  if  is  much  smaller  than 

"time  between  collisions",  the  change  in  v  will  be  almost  entirely  due  to  the 

o 

large  number  of  small  changes  due  to  distant  collisions,  the  probability  of  a  close 

collision  being  extremely  small.  The  net  result  of  a  large  number  of  small  effects 

generally  has  an  approximately  Gaussian  distribution  and  so  we  shall  assume  further 

a  Gaussian  distribution  for  v  .  That  is  we  assume  that  if  > 

o 
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gf(V0#t)  -  6(vo  -  Vq)  -  h  (Vq,0)  say 

Then  P  f(v  ,t  +~C  )  ■ 

°  0  '/ffTTrT 


(25) 


(v  -  v )  /f  m ....  ~ 

e  o  o  »h  (V  #XT  )  say  (26), 


where  the  "width" 


V7£7 


has  yet  to  be  determined.  In  fact  F('C)  will  depend  also 


on  but  since  we  are  only  considering  the  effect  of  collisions  on  a  localized 

distribution  (  A f(v  )  of  equatinn  (14)  we  can  ignore  the  dependance  of  F  on  V  . 

O  _  3 

We  now  note  that  a  physically  real  distribution  Of (v  »t)  can  be  decomposed 

'  o 

into  the  set  of  "modes"  h  defined  in  equation  (25)  i 


6f(vt)  ”  Sf(vt}  *  h{v0) 

From  equation  (26)  we  then  see  that 

-  F(v3|t)  »  h(v  fV  ) 


(27) 


(28) 


(Provided  we  neglect  collisions  between  pairs  of  the  electrons  "belonging"  to 
£)  which  are  negligibly  few  composed  with  those  between  dn  electron  from  S;f 
an  electron  of  4-  . 

Applying  this  result  twice  we  see  easily  that 


f(v  ,  t  +  2f)  -  f(vft+T)  »  h(v  ,  V  ) 

o  o  o 

■  f(v  ,t)  *  h(v  ,  "C  )  *  h(v  f  'C  ) 
oo  o 

a  f ( V  ,t)  #  h(v  ,2 ~C  ) 

0  o 

Whence  h(v  ,  2  'TT  )  ■  h(v  ,1P  )  •  h(v  ,TT  ) 

O  0  0 

The  right  hand  side  is  readily  calculable  and  we  obtain  finally 


F  (2  X  )  -  2  F  (  'tT  ) 
In  the  same  way  we  can  obtain 

F  (nt).nF  (rj 
and  so  F ( )  must  be  of  the  form 


f  ( 'tr )  -  2 
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where  is  a  constant  which  we  shall  leave  undetermined.  This  form  for  f*(  ) 

is  a  direct  consequence  of  the  assumption  of  the  Gaussian  distribution  for  0  v  . 

In  fact  06  is  just  the  of  Spitzer£T2]  ,  who  gives  formulae  for  its 

calculation.  Typical  values  in  the  eclar  corona  (n  ^  10^  cm~^,  T'-VlO  C,  V  5V ) 

15  2  -3  ®  so 

would  be  06^/3.10  cm  sec  .  Increasing  \IQ  by  a  factor  of  4  reduces  CL  by  almost 

2  orders  of  magnitude  (Indicating  that  the  assumed  Gaussian  ajipead  can  only  be  ♦••e'tfup 

for  very  nerrow  peaks). 

If  we  now  substitute  ^f(vQ)  from  equation  (14)  for  0  f (vQ,t)  in 
equation  (28)  the  result  is  i 


6" 


-i 


(v0  -  w  ) 


f(v  ,t  +r  )  -  -7==^==—^  e  2  y? 

0  V/27T  («?r+  w2  )  0 


(29) 


In  other  words  our  assumptions  are  such  that  if  0  f  is  initially  gaussian 
it  will  remain  gaussian,  but  it  will  spread  steadily.  Intuitively  we  know  that  this 
is  qualitatively  correct  except  that  the  mean  velocity  will  change  slowly  as  well. 

Also  if  the  amplitude  of  the  wave  is  sufficiently  large  that  the  number  of  trapped 
particles  is  large  the  rather  complicated  exchange  process  between  trapped  and  un¬ 
trapped  particles  will  provide  an  extra  source  of  instability  which  we  shall  avoid 
by  considering  only  small  amplitude  waves. 

With  the  reserves  in  mind  we  can  investigate  what  will  be  the  conse¬ 
quences  of  these  collisions.  The  equation  (2)  which  we  have  interpreted  as  a  dis¬ 
persion  relation  can  equally  well  be  interpreted  as  necessary  condition  on  f 

o 

for  thB  stability  ef  a  wave  of  given  k.  If  the  wave  is  stable,  k  will  not  change, 

and  if  initially  f  (v  ,0)  a  O/  f  +  f(v  ,0)  then  at  time  TT  later  we  have  seen 
°  o'  ^  eq  o' 

that  f  (v  )  m0(f  +Cf(v,t) 

□  o  eq  o  o 

2 

If  k  (tT)  calculated  from  equation  (2)  by  replacing  0f(vQ)  bV 

0  is  not  equal  to  the  original  k^,  (  k^(0)),  then  the  wave  is  unstable. 

This  will  almost  always  be  the  case.  However  Montgomery  [j has  ch*wn  that,  at 

least  in  some  cases,  the  solutions  of  B.G.K.  are  "asymptotically"  stable  -  and  so 

we  can  hope  that  if  k  (  "CT  )  k  (0)  then  the  waves  will  remain  fairly  stable. 

On  the  other  hand  thore  can  be  little  doubt  that  if  k^(^CT)  is  very  different  from 
2 

U  (U)  tlm  wave  will  be  very  unstable,  perhaps  evolving  towards  a  more  stable  state 
or  damping  out  completely. 

We  have  already  seen  that  if  the  maximum  of  0f  is  very  far  from  the 
wave  velocity  (more  precisely  if  w^  k^  )  the  width  of  the  gaussian 
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2 

curve  does  not  enter  the  approximate  expression  for  k  and  so  the  spreading 

of  the  gaussian  curve  will  probably  not  lead  to  instabilities;  in  this  case  the 

waves  approMimetely  satisfy  the  dispersion  relation  obtained  by  setting  f  «  f  . 

o  eq 

On  the  other  hand  when  w  and  W  are  both  email  (corresponding  to  a  sharp  peak  in 
the  velocity  distribution  near  the  wave  velocity)  k^  CZ,  )  (  calculated  on 
the  assumption  that  the  wave  is  stable)  will  chapge  very  rapidly,  leaving  little 
doubt  that  the  wave  is  unstable. 

For  example  if  w  *  0 


8  k2  (  T  ) 


%2 


n 

e 


€n 

w*  +  0(X 


(30) 


which  clearly  changes  very  rapidly  if  W2  /  OC 


is  small. 


These  unstable  waves  are  just  those  which  differ  greatly  from  the  dispersion 

relation  obtained  by  setting  f  =  f  . 

o  eq 

We  can  summarizei  and  at  the  same  time  generalize  these  results  as  follows. 

P  2  p 

If  does  not  imply  that  k  kEq  (  k^  that  given  by  setting 

fQ  ■  f  ).  However  the  difference  in  k2  is  only  due  to  very  high  first  (or  second, 

or  higher  ?)  derivatives  of  f  in  the  close  neighbourhood  of  the  wave  velocity. 

If  we  did  not  insist  that  f  be  continuous,  discontinuities  of  f  could  also 

o  _  o 

2  2 

introduce  large  differences  between  k  and  k  .  However  it  is  easy  to  see.  even 

eq 

without  making  the  above  calculations,  that  collisions  will  smooth  all  such  "roughness" 

2 

of  f  very  rapidly,  and  sincp  k  is  so  very  sensitive  to  these  details,  the  wave 
will  become  unstable  in  a  time  very  much  shorter  than  the  time  between  collisions. 

Hence  only  those  waves  for  which  f  is  very  "smooth"  in  the  region  of  the  wave 
velocity  can  conceivably  be  stable,  and  those  waves  approximately  satisfy  a  dispersion 
relation* 


Wb  shall  now  try  to  make  this  statement  a  little  more  precise,  at  least 
for  the  case  of  ^f  gaussian.  It  should  be  emphasized  however  that  what  follows 
depends  to  a  large  extent  on  the  additional  assumptions  made. 

For  instance,  to  make  more  precise  our  condition  that  fQ<^ f^  we 

shall  write  t 


|S f|  ■  b  f  (w  —  v  ) 

I  I  max  W  eq  p 


(31) 
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This  corresponds  to  the  case  where  the  initial  disturbance  effects  all  the  electrons 
more  or  less  equally,  i.e.  there  is  not  selective  acceleration  of  a  small  group  of 
electrons  (making  up  but  just  a  sort  of  "stirring  up"  of  the  distribution 

function  which  results  in  the  "roughness"  which  we  have  idealized  by  a  nerrow  gaussian 
peak.  Once  again,  the  validity  of  such  an  assumption  can  only  be  discussed  in 
terms  of  an  excitation  process. 

For  the  case  of  small  w,  which  is  the  only  case  we  need  consider,  we 
shall  rewrite  (31 )  t 

-£-0-<b  —lu.-ivv2 

W  N  \l 


where  b  is  some  constant,  less  than  1 . 

As  a  condition  for  stability  we  shall  arbitrarily  suppose  that  it  is 
necessary  that  t 


'd 


'dz 


(6  k2  (T  ))  <  B 


(32) 


8  k2  (  )  is  a  complicated  function  of  W,w,  'Z  etc.  However  it 


can  be  shown  that  the  average  value  of  this  function,  over  a  short  period  lOt*Z  ), 
for  small  values  of  w,  is  approximately  ; 


<-2-  s  x2<T  )> 


0>2 
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Hence  we  require  that  % 
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j£Ln—  / 

w4  ^ 
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0^2  a. 


Bearing  in  mind  ,that  t 

S<A 


max 


CQp 


oe 


(32') 


(17) 


we  can  readily{,V%y'’^tudying  figure  6,  that  (31’)  and  (32')  impose  an  upper  limit  of 

|6(k2)|  (The  "max"  refers  to  the  maximum  obtained  by  varying  w  holding  Cn 
and  W  fixed  i  here  we  are  interested  in  the  effect  of  varying  6n  and  W  as  well, 
subject  to  the  restraints  (31')  and  (32')).  In  figure  6  the  regions  of  the  (  £  n,  W) 
plane  which  do  not  satisfy  (31*)  or  (32*)  are  shaded  out.  The  lines  of  constant 


)8oo 


( 


&n 


/max  "^2 
that  above  a  critical  value  for 


constant  )  are  drawn  on  the  same  plane  and  it  will  be  seen 


)S(k2)/ 


m£x 


the  curves  are  contained  entirely  in  the 
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forbidden  (shaded)  region.  This  critical  value  is  of  course  our  upper  limit  and  so  we  have 

|6(k2)|  A  exp  (-1/3  V2  /v2) 

[max  2  p 

where  A  is  a  parameter  depending  on  60  and  V  and  on  the  assumed  values  of  b 

(  J6fj  ^  bf  )  and  B  (  JiL-  0  (k2)  <  B). 

1  eq 

2  2  12' 

In  figure  (4)  the  curves  k  and  k  +  l£sk  j  are  sketched  as  functions 

“  'max 

of  v^  which  is  almost  exactly  equal  to  the  wave  velocity.  In  figure  (5)  these  curves  are 
re-drawn  in  terms  of  the  coordinate  used  by  Denisse  and  Delcroix  when  discussing  dispeiw 
sion  relations.  We  see  that  a  dispersion  relation  exists  for  those  waves  which  travel  much 
faster  than  the  thermal  velocity  and,  as  we  have  already  shown  jn  section  IV  this  disper¬ 
sion  relation  is  well  represented  by  approximate  dispersion  relation  of  Bohm  and  Gross 


* 


mik2 


It  should  be  noted  however  that  this  dispersion  relation  now  only  describes 
those  waves  which  are  relatively  stable  in  the  presence  of  collisions.  On  the  other 
hand  there  appears  to  be  no  relation  between  QJ  and  v  »  CO/k  when  v^  is  of  the 
order  of  V  the  thermal  velocity.  In  this  respect  it  should  be  noted  that  in  fact 
Denisse' s  equations  only  describe  the  potential  in  the  region  between  two  consecutive 
potential  minima.  To  describe  a  periodic  wave  it  is  necessary  to  assume  that  the  trapped 
electrons  have  the  same  distribution  between  each  pair  of  potential  minima,  an  assumption 
which  is  of  little  consequence  in  the  case  of  waves  travelling  faster  than  V,  but  which 
needs  to  be  reconsidered  for  waves  travailing  slower  than  V  (the  thermal  velocity).  In 
fact  it  is  not  hard  to  see  that  in  a  given  stable  wave  travelling  slower  than  V,  the 
distance  between  potential  minima  may  vary  from  one  minimum  to  the  next  and  the  potential 
observed  at  a  point  fixed  in  the  plasma  will  vary  more  or  less  at  random.  There  is  there¬ 
fore  no  real  reason  to  distinguish  between  such  waves  and  thermal  fluctuations  -  the  sta¬ 
ble  waves  of  low  velocity  being  like  a  special  case  of  thermal  fluctuations. 

All  the  above  calculations  have  been  made  for  the  very  special  case  of  a  dis¬ 
tribution  made  up  of  two  Gaussian  functions.  The  calculations  can  be  immediately  extended 
to  distribution  functions  made  up  of  a  finite  number  of  gaussian  functions.  Ii  is  possible 
that  the  calculations  can  be  extended  to  the  case  of  a  general  distribution  function.  This 
has  not  as  yet  been  attempted  as  it  is  doubtful  if  the  model  used  for  the  collision 


effects  justifies  an  elaborate  mathematical  theory.  At  the  same  time  the  physical  picture 

is  clear  t  for  those  small  perturbations  of  a  plasma  which  propagate  as  stable  waves  in 

the  absences  of  collisions  we  can  suppose  that  f  (v  )  s\sf  (v  -  v  ). 

o  o  eq  o  P 
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.  This  however  is  not  a  strong  enough  con¬ 
dition  to  determine  J<  .  However  those  f  (v^)  which  give  a  J<  greatly  different 
from  that  given  by  the  case  f  (v  ) 

*9*0  °  ° 


o  o 

f  (v  -v  )  are  characterized  by  large  values 
eq  o  p 


of 


in  the  neighbourhood  of  vq  ■  0. 


If  we  take  the  characteristic  effect  of  collision  to  be  the  introduction 
of  a  spread  in  the  velocities  we  see  that  such  critically  steep  slopes  cannot  exist 
for  any  length  of  time.  On  the  other  hand,  if  the  collisionless  theory  has  any 
meaning  at  all,  it  is  in  the  case  where  the  distribution  function  is  smooth,  and 
so  we  come  back  to  the  dispersion  relation  found  in  linear  theories  -  at  least 
when  the  wave  velocity  is  large  -  but  without  the  damping. 

The  points  are  i  i  Collision  damping,  or  rather  collision  induced 
instability  may  occur  in  a  time  much  shorter  than  normally  assumed. 

2  The  most  stable  waves  are  probably  those  which 
satisfy  the  linear  dispersion  relation. 
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SECTION  VII  -  CONCLUSION  AND  DISCUSSION 

This  report  is  largely  concerned  with  the  nature  of  stable,  small 
amplitude,  longitudinal  waves  in  a  collision  free  electron  plasma.  However,  since 
the  results  of  Landau  and  Van  Kampen  are  generally  taken  as  meaning  that  no  such 
waves  can  be  stable,  a  considerable  part  of  the  report  has  been  taken  up  by  argu¬ 
ments  concerning  the  existence  of  non-dampod  waves. 

After  this  the  approximate  results,  obtained  by  Denisse,  are  used  to 
examine  the  existence  of  a  dispersion  relation  for  these  stable  waves,  subject 
to  the  assumption  that  the  distribution  function  in  the  wave  is  everywhere  appro¬ 
ximately  equal  to  the  equilibrium  distribution.  Such  an  assumption  excludes  the 
interesting  case  of  a  wave  excited  by  a  beam  of  particles  injected  into  the  plasma 
but  since  we  have  as  yet  no  description  of  possible  excitation  processes  this 
appears  to  be  a  necessary  exclusion.  Even  subject  to  all  these  restrictions 
(stable  waves,  f  f  )  it  turns  out  that  there  is  no  dispersion  relation,  i.e. 

’  eq 

any  wavelength  can  correspond  to  any  wave  velocity  (at  least  for  wave  velocities 
less  than  the  velocity  of  light). 

However  in  Section  VI,  by  a  rather  intuitive  considerations,  it  has 
been  shown  that  the  effect  of  collisions  is  probably  to  render  highly  unstable 
those  waves  which  do  not  satisfy  approximately  the  normally  obtained  dispersion 
relation. 

The  work  described  is  of  course  strictly  limited,  for  example  the  general 
problem  of  the  evolution  of  a  plasma  after  an  initial  time,  attempted  by  Landau, 
has  not  been  attempted,  Nor  is  .'t  clear  how  it  can  be  approached.  It  has  been 
argued  by  Bernstein,  Greene  and  Kruskal  and  by  Weenink  that  the  B  G  K  "modes" 
cen  be  identified  with  Van  Kampen  modes.  However  the  Van  Kampen  modes  are  inde¬ 
pendent  in  that  they  satisfy  a  linear  differential  equation.  On  thB  other  hand 
there  is  every  reason  to  suppose  that  the  stable  waves  which  can  exist  in  a,  plasma 
cannot  in  general  be  considered  to  propagate  independently. 

Bohm  and  Gross  have  already  pointed  out  that  it  is  a  good  approximation 
to  consider  two  waves  propagating  in  opposite  directions  (or  more  generally  with 
highly  different  velocities)  to  be  independent.  On  the  other  hand  the  implication 
is  that  waves  with  almost  the  same  velocity  will  interact  strongly.  That  this  is  so 
is  easily  sesn  t  the  properties  of  the  wave,  such  as  the  wavelength,  depend  strongly 
on  the  distribution  of  electrons  near  the  wave  velocity.  On  the  other  hand  it  is 
just  these  electrons  which  will  be  most  strongly  perturbed  by  another  wave  with 
almost  the  same  velocity.  Hence  the  two  waves  certainly  interact  through  those 
electrons  close  to  both  wave  velocities. 


26 

The  identification  of  the  B  G  K  equilibria  with  the  Van  Kampen  modes 
does  not  therefore  automatically  justify  the  application  of  Van  Kampen'a  analysis 
for  which  the  independence  of  these  modes  is  essential. 

It  is  possible  therefore  that  the  analysis  of  stationary  waves  (including 
that  describes  here)  will  not  be  directly  applicable  to  a  study  of  the  more  general 
problem.  On  the  other  hand  the  concepts  exposed  here  and  elsewhere  should  be  of 
some  value  in  the  development  of  the  more  general  theory. 
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Poissons  equation  becomes 
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If  we  had  chosen//A<C’  0  this  would  have  changed  the  sign  of  both  terms  in£^"  "j ^ 
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and  xq  is  a  constant  of  integration. 


We  can  readily  plot  (x  -  x^)  as  a  function  of  z,  or  by  distorting  the  scale, 
In  other  worts  we  have  x  as  a  function  of  Cp  the  inverse  function  of  Cp  (x). 

We  can  easily  see  that  Cj 9  is  an  analytic  function  of  x.  It  is  true  that  we  have 
ignored  the  motion  of  the  ions,  but  there  is  no  reason  to  suppose  that  this  intro¬ 
duces  any  fundamental  difference. 
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We  wish  to  study  the  function  t 
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By  definition  the  Hilbert  transform  of  the  function  f(x)  is  proportional 
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We  shall  therefore  make  use  of  the  general  rules  concerning  this  transform, 
which  (except  for  the  first)  can  readily  be  justified  from  the  definition 
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It  can  readily  be  verified  that  the  general  solution  is 
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which  is  the  required  result  (equation  (12)) 
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x  e  dx  «  0  if  r  is  odd 


rco 


r-1 

t  2  dt  if  r  is  even. 


th 


Since  the  n  term  of  the  series  in  equation  (Cl  )  is  greater  than  the  (n  -  1) 
2 

term  if  n  >  y  +  ^  it  is  clear  that  the  series  diverges,  and  therefore  that 


th 


fL  (y)  - >  +  OO 

m  — >co 


We  intend  to  show  that  the  expansion  of  equation  (Cl)  is  nevertheless 
useful  as  on  asymptotic  series. 


xxxxxxxxxxxxxx 
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If  we  write 


and 


m  i  \  n  -x 

f  lx)  ■  x  e 
n 


then 


R2{m-1)  ”  5T3"  9 2m  W 

y 


Proceeding  exactly  as  in  appendix  B  it  is  easy  to  show  that 


(C  3) 

(C  4) 


g2m(y)  ■  2P  (m  +i)  y2"1  e"y 


u  ~2m  .  r  v 

e  u  du  +  C  r 


(r  5) 


where  _a_  ->  0  (since  we  are  interested  in  the  case,  y  0)  and  in  principle  C 

can  be  determined  for  given  a.  For  example  it  can  be  seen  intuitively  that  g  (y) 

“  2m 

■  0  for  some  y  ■  y_  ^-v.  \frni  we  could  write  a  ■  y_  and  C  ■  0. 

■m  2m 

It  suffices  however,  to  note  that  ^  j  *  ~  ~  (_^<  +  OO  monotonically  and 

is  therefore  always  positive  for  some  y  >  Y  (m)  which  we  need  not  evaluate. 

Hence  R^fy)  0  for  y  Y(m)  (C  6) 

But  from  (C  1 ) 


i.e. 


Vy)  “  R2m-2(y)  -  + 


R2m(y)  ^  R2m-2  (y) 


r<~ 


(m  +  \) 


2m 


>° 


(C  7) 


Combining  (C  6)  and  (C  7)  we  see  that  for  y>Y(m) 
R2m(y)!  <[R2m-2(y)  / 


If  we  can  invert  the  relation  y  ^Y(m)  to  read 
m  ^M(y) 


(C  0) 


then  M  (y)  gives  us  the  number  of  terms  to  take  for  the  best  approximation  to 

91 (y)» 

If  the  intuitive  result  mentioned  above 
Y  (m)  ryj  y/nT11 
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ia  correct,  then  y  Y  (m)  m<  y  (approx.) 


i«B* 


M  (y)  O—  y 


(which  corresponds  approximately  to  the  smallest  term  of  the  series). 

It  can  be  shown  that  R_  (y)  . . 0 

2m  . 

y  — >OQ 

Choose  a  ^y/’m'  so  that  eU^  ui  eu2 

for  all  Uj  ^  a 

-y2  i  /”  y  u2  -2m 


Then  y2™  e  ^  e"  u  "'  du  +  C 

2m  -y2  2m  -y2  ["  ^  u2  -2m 

-ye  C  +  y  e  J  e  u  du 


J 


2m  -y' 
+  y  e 


2  fy 


2  -2m  J 

u  u  du 


Ay 


The  first  term  tends  to  zero  when  y-^QQ(0  <  A<^1) 
When  y  a/A  the  second  term  is  I 


2m  -y 


2  r  Ay  2 


/ 


2  /-Ay  (».,>2  v-2m 

J  0 


u~  -2m  ^  2m  -y“  f  ”J'  (Ay)  (Ay)  du 

e  u  du  ^  y  e  «  «  '  * 


(Ay  -  a)  e-(1  -  A2)  y2 


2m 


The  third  term  is 
2m  -y2 

y  e  J 


Ay 


A 

— ^  0 


2  __  2m 

u  -2m  ^  y 

e  u  du  ^ 


{  y->oo  ) 

_y 2  r  y 


e  u  du 


Ay 


-y 


TrTTT 


FJ 


y 

— }  o 

Ay 


(y  — >CO) 


2m. 


Clearly  the  R^ly)  0  ^R2m/"'~'1^y  )  (y 


2m 


The  functions  R__(y)  are  sketched  in  figure  (Cl). 
2m 
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APPENDIX  D 


Landau's  result  is  that  asymptotic  behaviour  in  time  of  a  plasma  is 
characterized  by  a  potential  wave 

i(kx-O)t)  (D  1) 


where  k(real)  and  QJ  ( complex )  satisfy 


CO 


F»  (v) 
o 


v  -  aYk 


dv 


(D  2) 


where  the  path  of  integration,  C,  is  that  indicated  in  figure  (D  1),  which  passes 

below  the  pole  (v  ■  CO/k)  of  the  integrand.  The  function  F  is  supposed  to  have 
o  o 

been  continued  analytically  into  the  whole  complex  velocity  plane.  We  have  used  F^_ 

in  (D  2)  rather  than  Landau's  f  ,  to  distinguish  from  Denisse's  f0 

We  shall  show  here  that  we  can  choose  F  to  satisfy  equation  (24),  such 

o 

that  there  exist  real  CO  and  k  satisfying  equation  (D  2). 

Suppose  that  CO»k  (both  real)  and  F  satisfy  equation  (D  2).  In  this 


case  we  can  replace  Jc  by  Jr 


r 


L 


where 


f~  ,  Y  are  the  contours 


shown  in  figure  (D  2) 

The  assumption  that  F'  can  be  analytically  continued  into  the  complex 

°  i©  V 

plane  means  that  on  the  semi-circle  v  «*  CO/  k  +  r  e  which  the  contour  o 

definies  we  have 


F*  (  00/k  +  r  e1^  )  -  F»  (  CO/k)  +  r  e1®  P*  (  CcVk)  +  0  (r2) 
o  o  o 


v  =  CO/k  +  r  e ^  -  i  r  e*® 


d  e 


ILM 


FA  (CO/k+reiQ) 


Fj  (QUA) 


v  -  on 

CA>k  + 

f  :l.  ±l 

r  s 


iQ 


/VC  I 


v  -  CO/k 


r_njoo/k) 


+  F"  (  CU  /k)  +  0  (r) 
o 


+  F"  (CU/k)  +  0  (r) 
o 


A p.  0  2 


-  i  F  (  CxJ/k) 
0 


r2.1T 

I  d  9  +  0(r) 


i77>  tOJ/k)  +  DM 


Hence 


I  ~ 

J S*  V  - 


(v) 

dv  - 

OO/ k*  J 


(tfk  -r  /-oo 

Fn  (V>  • 

- 2 -  dv+  i  F  (C*A)  +  0(r) 


+ 

— do 


V  -  COA 


Since  the  left  hand  side  is  independent  of  r,  and  the  equation  is  true  for  any 
r  0,  we  may  take  the  limit  ae  r— >+0. 


F  (v) 
o 


'°°  F*  (v)  , 

dv  -  P  I  — -  d*+  i(Fo(v))v  -COA 


(D  3) 


,  v  -  OJA  J  v  -  COA 

'C  s  J~O0 

This  equation  is  only  exact  if  COA  is  real. 

2  * 

We  now  see  the  condition  that  k  be  real  requires  that  (F  (v))  ,tVl 

o  v*«  tif  k 

be  zero. 

The  problem  is  therefore  to  choose  d  n  F  (v) ,  which  can  be  continued 

°  * 

analytically  into  the  whole  complex  plane,  which  satisfies  (v)  ■  0  for  some 
v  and  which  satisfies  relation  (24),  i.e. 


F  (v)  -  f  (v)  I  <C  $L 
a  eq  I 


(24) 


These  conditions  can  clearly  be  satisfied  by  writing  s 


F  (v)  -  f  (v)  +  6f<v) 

o  eq 

.  f  v  - v  \2 

where  f(v)  *  OCe  V  /$  / 

subject  to  the  condition  that  j  0(.  j  ^ 

Then  ^  can  always  be  chosen  small  enough  that  F^  (v)  will  have  a  zero 
for  some  v  ■  v  . 


To  ensure  that  f  F  (v)  dv  ■  If  (v)  dv  we  must  replace  D  f(v) 
-  Jo  J  eq 


by 


i  £  F(v)  -  i  S>  f  (2v0  -  v)  which  is  an  odd  funtion  of  v  -  v  .  Clearly  this  does 


not  change  the  derivative  s  F  (v  ). 

oo 
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This  suffices  to  prove  that  there  exist  a  real  CO  and  a  real  k  which  satis¬ 
fy  equation  (D  2),  subject  to  the  relation  (24),  for  any  ^^>0. 


To  (v) 
J  V  -  V 


dv  {  CO  «*  vQ  k) . 


STABLE  WAVES 


FIG  4*  FIG 
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